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CONSTANT  RADIUS  MAGNETIC  ACCELERATION  OF  A 
STRONG  NONNEUTRAL  PROTON  RING 

I.  Introduction 

As  a result  of  an  intense  experimental  effort  during  the  last  year, 
pulsed  ion  beams  are  presently  available1?2  at  a peak  power  level  in 
excess  of  0.2  Terawatts  at  energies  of  about  1 MeV.  However,  several  of 
the  fusion  related  potential  applications  of  ion  beams,  including  target 
irradiation  experiments3  5 and  field  reversing  ion  ring  fusion  reactors6 
require  power  levels  in  excess  of  100  Terawatts  at  energies  10-10,000 
MeV.  Thus,  it  is  desirable  to  develop  methods  of  acceleration  that 
would  eventually  increase  the  energy  of  these  intense  ion  beams  by  as 
many  as  four  orders  of  magnitude. 

Free  or  adiabatic  magnetic  compression  of  electron  rings  at  the 
University  of  Maryland7  led  to  a two  hundred  fold  increases  of  the 
electron  energy.  The  application  of  this  acceleration  method  to  ion 

o fl 

rings  has  been  considered  by  Fleischmann  and  Sudan  and  Ott  . In  the 
adiabatic  compression  the  magnetic  field  varies  with  time  and  thus  the 
magnetic  moment  (m.  = Mqv^/2B,  nonrelativistic)  of  a gyrating  charged 
particle  is  an  adiabatic  invariant.  Therefore,  the  energy  of  the  com- 
pressed particles  increases  linearly  with  magnetic  field.  The  final 
value  of  magnetic  field  required  to  accelerate  a 1 MeV  proton  ring  of  1 
meter  initial  major  radius  to  about  1 GeV  is  in  excess  of  1.4  M Gauss. 

In  addition,  extraction  and  unwinding  of  the  ring  after  compression  into 
Note:  Manuscript  submitted  March  10,  1977. 
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a low  divergence  straight  beam  is  rather  difficult.  The  conservation  of 
the  canonical  angular  moment™  Pg  warranties  a rapid  radial  expansion  of 
the  beam  when  a ring  made  up  of  ions  with  Pg/^o  is  extracted  axiaJJy 
through  a half-cusp.  Radial  extraction  of  the  beam  is  also  made  diffi- 
cult because  the  reduction  of  the  beam  major  radius  with  increasing 

magnetic  field  excludes  the  use  of  deflecting  electrodes. 

By  contrast  to  the  free  compression,  in  the  constrained  or  constant 

radius  magnetic  acceleration  the  magnetic  moment  is  not  a constant  of 
the  motion  and  the  energy  of  a nonrelativistic  gyrating  charge  particle 
increases  proportionally  to  B2.  In  addition,  radial  extraction  of  the 
beam  can  be  accomplished  using  magnetic  or  electric  deflecting  fields. 

The  radius  of  the  ring  may  be  kept  constant  either  by  the  applica- 
tion of  a positive,  radial  electric  field  that  increases  with  time  or  by 
shaping  the-  external  magnetic  field  in  such  a way  that  the  betatron  flux 
condition,1'1’12  i.e.,  < B (t)  > = 2Bo(t),  where  < B(t)  > is  the  average 
and  BQ(t)  the  magnetic  field  at  the  orbit  of  the  particle,  is  satisfied. 

In  the  first  approach,  the  value  of  the  radial  electric  field  re- 
quired to  keep  the  radius  of  ions  constant  depends  u^£  the  initial  and 
final  value  of  the  axial  magnetic  field  and  the  initial  azimuthal  veloc- 
ity of  the  ions.  The  value  of  the  maximum  field  required  for  the  accel- 
eration of  a 5 m radius  ring  from  2 to  1,000  MeV  is  about  2.5  MeV/cm. 

This  field  is  by  a factor  of  ^ higher  than  the  vacuum  breakdown  field. 
However,  the  breakdown  in  the  vacuum  gap  may  be  avoided  because  the  ex- 
ternally applied  axial  field  is  higher  than  the  critical  magnetic  field13 
required  for  magnetic  insulation.  Details  about  this  method  of  accelera- 
tion will  be  reported  in  a future  publication. 
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In  this  report  we  present  theoretical  results  on  the  constant 
radius  acceleration  of  a strong,  nonneutral  ion  ring  in  a modified  beta- 
tron field.  Superimposed  on  the  axial  magnetic  field  that  satisfies 
the  betatron  flux  condition  is  a strong  Bg  field.  This  magnetic  field 
configuration  was  initially  used  in  the  plasma  betatron. 14-16  However, 
the  orbit  stability  analysis  in  the  plasma  betatron  has  been  restricted14 
to  weak  electron  beams.  As  it  is  shown  later  on,  the  Bg  field  is 
necessary  for  the  stability  of  the  orbits  when  the  self  fields  of  the 
ring  are  not  negligible. 

At  power  levels  ~ 1012  watts,  the  constant  radius  induction  ion 
acceleration  competes  with  several  collective  ion  acceleration  schemes. 
One  such  scheme  is  the  ionization  front17  19 acceleration  process  in 
which  an  intense  electron  beam  is  used  to  form  a single  accelerating 
electrostatic  potential  well  which  traps  and  accelerates  a group  of  ions. 
By  carefully  controlling  the  ionization  rate  the  ions  can  be  accelerated 
in  a controlled  manner  to  energies  many  times  the  electron  beam  energy. 
Another  class  of  accelerators  rely  on  the  controlled  acceleration  of 
ions  trapped  in  a wave  which  is  supported  by  an  intense  relativistic 
electron  beam.  The  Auto  Resonant  Accelerator^0  (ARA)  is  of  this  type. 

In  the  ARA  the  ions  are  trapped  in  a beam  cyclotron  wave  which  is 
accelerated  on  an  intense  relativistic  electron  beam  by  gradually  re- 
ducing the  guiding  axial  magnetic  field.  In  the  Converging  Guide 
Accelerator21  (CGA)  a space  charge  wave  on  an  intense  relativistic 
electron  beam  propagating  through  a cyclindrical  waveguide  is  accelera- 
ted by  adiabatically  decreasing  the  waveguide  radius.  In  both  the  ARA 
and  CGA  scheme  the  final  phase  velocity  of  the  accelerated  wave  is 
roughly  equal  to  the  electron  beam  velocity. 
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II.  Induction  Accelerations 

A low  density  monoenergetic  ion  ring  in  a locally  uniform  external 
axial  magnetic  field  will -execute  gyro-motion -around- its -guiding-center-. 
The  particle  radius  is  given  by 


:o  - ^(t)/  [no(t)/vo0(t>]  , 


(1) 


where  v Q is  the  particle  velocity  perpendicular  to  the  slowly  changing 


local  magnetic  field  B , |e|Bo/Moc  is  the  nonrelativistic  ion  cyclo- 

1 

tron  frequency  and  r , „-i-5 

Vo6  - HV')2]  • 

The  aximuthal  electric  field  at  radius  r,  induced  by  the  changing 
axial  magnetic  field  is 


§£  "ft 


(2) 


where  (b>  is  the  average  value  of  the  magnetic  field  within  the  radius  r. 
The  rate  of  change  of  the  total  energy  W of  a charged  particle  of  mass  Mq 
and  charge  q = je|  in  the  presence  of  an  electric  field  E0  is  given  by 


dW  | i ti 

— = I e]  v.E  =|e|voQ  E0, 


(3) 


where  W = Mq  c2Yqq . Substituting  Eqs.  (l)  and  (2)  into  Eq.  (3),  it  is 


obtained 


(»0 
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Equation  (l)  is  valid  only  for  slowly  changing  fields,  i.e.,  when 


[in  (Bq)] 


< < n . 

o 


Since  Eq.  (l)  has  been  used  in  the  derivation  of  Eq.  (4),  this  equation 
is  also  valid  when  the  inequality  in  Eq.  (5)  is  satisfied. 

In  the  case  of  a uniform,  externally  magnetic  field  i.e., 

Bq  (t)  = < B(t)  >,  Eq.  (4)  gives 


Y a (o)  +1  B (t) 

'!<*>  ■ kfwH  kror  ?i(o)’ 


where  £\(o)  = ^Yqq (o)  - lj  Mqc2  is  the  initial  ion  energy,  £\(t)  = 

|^Yoe ( t)  - lj  Mqc2  is  the  final  ion  energy  and  Bo(t)/BQ(o)  is  the  ratio 
of  the  final  to  the  initial  value  of  magnetic  field.  -Combining  Eqs.  (l) 
and  (6),  the  major  radius  of  the  ring  r (t)  is  given  by 


■o(t)  = ro(o)  [Bo(°)/Bo(t)] 


Therefore,  in  the  free  or  adiabatic  acceleration  the  energy  of  the  non- 

relativistic  ions  increases  proportionally  to  the  applied  magnetic  field 

i 

and  the  major  radius  of  the  ring  varies  as  Bq2. 

By  contrast,  when  the  average  magnetic  field  is  twice  the  local 
field,  i.e.,  when  < B(t)  > = 2B  (t),  Eq.  (4)  gives 


" V e (°)  + 1 1 B (t)  "I  2 

,<t)=  [v°e(t)  + i J <iw 
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Hence,  the  ion  energy  in  the  nonrelativistic  case,  increases  as  the 
square  of  the  magnetic  field.  This  is  a significant  advantage  over  the 
free  compression  scheme  where  the  energy  [Eq.  (6)]  is  proportional  .to  the 
magnetic  field.  In  addition,  Eq.  (8)  shows  that  the  ratio 
c2  |\20(t)  -1  J /^2(t)  remains  constant  during  compression.  However, 
this  ratio  is  according  to  Eq.  (l)  equal  to  the  square  of  the  ring 
radius.  Therefore,  when  the  betatron  flux  condition  BQ(t)  = < B(t)  >/2 
is  satisfied  the  radius  of  the  ring  remains  fixed,  although  the  magnetic 
field  increases  with  time. 

In  the  conventional  betatron  accelerators  the  required  spatial 
variation  of  the  magnetic  field  is  achieved  by  introducing  ferromagnetic 
material  into  the  core  region  of  the  system.  In  a proton  accelerator 
the  presence  of  ferromagnetic  material  is  not  desirable.  A ferromagnetic 
core  in  a proton  accelerator  would  present  rise  time  and  saturation  limi- 
tations on  the  field  as  well  as  weight  problems.  However,  a field  con- 
figuration having  the  proper  radial  nonuniformity  can  be  easily  be 
achieved  with  air  coils. 

III.  Orbit  Stability  in  a Modified  Betatron  Field 

In  this  section  we  analyze  the  stability  of  the  single  particle 
orbit  in  a modified  betatron  field.  In  such  a field  configuration,  an 
aximuthal  time  independent  magnetic  field  B0  is  superimposed  on  a time 
varying  betatron  field,  such  that,  at  r = rQ  the  condition  Bo  = i < B > 
is  satisfied.  The  present  analysis  is  not  restricted  to  weak  rings  and 
thus  both  the  self  electric  and  the  self  magnetic  fields  are  included. 

For  rings  with  appreciable  self  fields,  stable  orbits  exist  only  if  the 
B0  field  is  considerably  higher  than  the  initial  betatron  field  B . 
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As  shown  in  Fig.  1,  the  Bg  field  is  conveniently  produced  by  passing  an 
axial  current  along  the  symmetry  axis  of  the  system  and  the  betatron 
field  is  generated  by  two  concentric  air  coils  having  different  diameters. 
For  the  analysis,  the  external  fields  are  taken  as 


Bz(r,t)  = Bo(t)(ro/r)n  , 


B ( r , z , t ) = -nB  (r,t)z/r  , 

L 


Bfi (r)  = B r /r 
W N 1 o ? 


1 r &B  (r',t)  , 

SQ(r,t)  = - — / ^ r ^ 


(9a-d) 


where  n is  the  external  field  index,  Bg(r)  is  the  field  produced  by  the 

axial  current  I , EQ(r,t)  is  the  induced  electric  field  and  (r  , z ) 
o 7 u o o 

are  the  coordinates  of  the  average  position  of  the  particles  forming  the  _ 
ring.  Note  that  Eqs.  (9a, b)  are  the  expressions  for  the  local  field  and 
n is  assumed  constant  and  uniform  over  the  ring. 

To  approximate  the  effects  of the  self  fields  we  replace  the  ring 
by  a cylinder  of  radius  a and  length  2rrrQ.  The  magnitude  of  the  total 
force  due  to  the  self  electric  and  magnetic  field  a distance  £ from  the 
axis  of  the  cylinder  is 


F (t)  = W nj(t)  p 

S(t)  2 Y^TtT  5 ’ 


where  m(t)  is  the  ion  density  in  the  ring,  which  is  assumed  uniform. 
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The  total  self  fields  of  the  ring  can  be  represented  by  an  equivalent 
self  electric  field  of  the  form 


-*  4n|e|n  (t)  r * , a 

"self^2’^  = 2V0?'(tT  L(r"ro)  ^r  + C«o>  ez  J 


where  e and  e are  unit  vectors  in  the  r and  z directions.  Equation 
r z 

(ll)  is  a good  approximation  for  large  aspect  ratio  rings,  i.e.,  when 
(r0/a)>  > -1" 

We  now  perform  a perturbation  analysis  using  the  external  fields  of 
Eqs.  (9)  and  the  equivalent  self  field  of  Eq.  (ll).  The  particle  veloc- 
ity components  are  written  in  the  form 


v (t)  = v + 6v  , 
r or  r 


l(t)  = v Q(t)  + 6v0 


v (t)  = v + 6v 
zv  ' - oz  z 


where  v = v =0  and  the  perturbations  &v  , 6vfi  and  6v  are  much 
or  oz  r ° z 

smaller  than  v0.  Substituting  the  above  velocity  components  into  the 
relativistic  orbit  equations  and  keeping  terms  to  first  order  in  the 
small  perturbations  we  obtain 


(Y„fi6v  ) + Yr 


dt  ' To0  r 


C-'f 

Vw 


(l-n-n  ) 6r 
' s 


= -^1  6vz  ■ no  YO0  6ve 
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(y3  6Vfl ) = ) 

dt  'Yo0  & M ' dr  J 


&r 


r = r 


— (y  a 6v  ) + Y o ( — — ' ^ (n-n  )6z 
dt  woO  z'  T o6  \ Yo0  / s' 


= 0i  5v  , 
r 


(I3a-c) 


where  Y0q  - 


[l-(vo8(t)/=f]  -4 
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"s  f2Y8(t)  no(t>] 


is  defined  as  the  self  field  index. 


ou  = (4tt  e2  ni(t)/MQ)‘s  ,^i  = |e|Bi/MQc)  and  6 r,  6z  are  small  spatial 

perturbations  about  rQ,  zq,  i.e.,  d6r/dt  = Sv^,.  It  can  be  easily  shown 

from  Eq.  (9d)  that  (dEg/dr)  = o when  the  betatron  field  condition  is 

r - r 


o 

satisfied.  Hence  from  Eq.  (13b)  we  see  that  the  perturbed  angular 
velocity  6v0  vanishes.  Equations  (13a)  and  (l3c)  can  now  be  put  into 
the  rather  simple  form 


6r  + 


Yo0 


6f  — uu2  6r=- 
r 


6z 


ro9 


(l4a,b) 


and 


+ -2^-  6z  + u)2  6z  = 


0, 


Y 


o6 


Yo0 


> 


where  a dot  represents  the  operator  d/dt,  uj2  = (^0/voq)2  (l-n_ns)  an(i 

rj2  _ Cn  /y  * )2  (n-n  ).  The  above  equations  for  the  small  displacements 
z “ o ow  s 
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r 


6r  and  6z  about  r and  z are  coupled  through  the  azimuthal  magnetic 
o o 

field.  The  quantities  Yqq , and  ns  are  assumed  to  be  "weak"  functions 

of  time  compared  to  the  cyclotron  period. 

The  self  field  index  n , can  be  expressed  in  terms  of  the  field 

s 

reversal  factor  £ = 1 Bself/Bo^  where  Bself  and  Bo  iS  the  Self  and  ex“ 
ternal  magnetic  field  at  the  inner  surface  of  the  ring.  That  is, 

^pi  ro  C , (15) 

113  ' = a Y20f32e 

where  P Q = v Q/c.  Although  £ may  be  small  (a  few  $ field  reversal)  the 
o»  o» 

self  field  index  can  be  large,  ns  » 1.  Since  the  magnitude  of  the 
external  field  index  is  typically  of  the  order  of  unity,  [ nf  = 0(l),  the 
radial  and  axial  betatron  frequencies  are  approximately  equal,  that  is, 

w2  = m2  s uj2 (t)  < o.  (16) 

r r o' 

Using  (l6),  the  condition  for  the  orbits  to  be  stable  is  shown  in 
appendix  A to  be  given  by 

e2(t)  < 1,  (17) 

where  e2(t)  = -^/(C^Ay  q)2  > o.  The  stability  condition  in  (17) 
written  in  terms  of  the  self  field  index  is  ng  <(n1/^Q)2A. 

In  the  absence  of  a B0  field  the  orbits  are  stable  if  o < n < 1 

and  n < h.  The  total  number  of  ions  and  current  in  the  ring  is 

s 

given  by 
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N.  = 2 X 1016  (— ) v3  «p2e  n 
ion  V o»  o#  s 


I.„„  = 1.45  X 1 07  (^-)2 


0eF  O0  s 


(l8a,b) 


where  a and  r are  in  centimeters  and  I.  is  in  amperes.  Since  the 
o ion 

field  reversal  index,  number  of  ions  and  current  in  the  ion  ring  is 
proportional  to  the  self  field  index,  it  is  desirable,  for  stability- 
reasons,  to  have  initially  as  large  as  possible. 

We  now  proceed  to  analyze  the  solution  of  the  linearized  coupled 
orbit  equations  Eqs.  |(l4  a,b)J  when  ng  > > 1 and  | n|  = 0(l).  To 
facilitate  the  analysis  we  further  assume  that  e2  < < 1.  The  solutions 
of  Eq.  (l4)  under  the  above  assumptions  (see  appendix  A for  a general 
derivation)  are 


J 6(t')dt'V-  _i  / (o)/. 's  ' 

>r  = e V J n (t  )dt 


6z  = e 


/ n O 


(I9a,b) 


where 


6(t)=  ife(e2) 


> (20) 


a , b are  constants  determined  by  the  initial  conditions  and  u>  ' 
n’  n n 
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I 

I 


the  roots  of  the  zeroth  order  dispersion  relation 


h(o)> 


2 2 

- 


(t)  n'- 


Since  the  frequencies  m K J are  all  real  the  radial  and  axial  dimensions 

^ t "1 

of  the  ring  vary  as  exp  £ 6(t  )dt  « 1 + |^e2(t)  - e (°)J  /2. 

The  cross  sectional  area  of  the  ring  at  time  t is,  therefore,  given  by 


A (t)  * A(o)  {l  + e2(o)  [e2(t)/e2(o)-l]}  , (22) 

where  A(o)  is  the  initial  cross  sectional  area.  Using  the  definition  of 
e2(t)  and  the  fact  that  the  total  number  of  particles  in  the  ring  is 
conserved  we  find  that 


!(t)/e2(o)  = [a(o)  YoO(o)/A(t)Yo0(t)] 


Combining  (22)  and  (23),  the  cross  sectional  area  as  a function  of  time 


is  given  by 


Thus,  the  area  A(t)  decreases  as  the  magnetic  field  increases  and  is 
bounded  by  the  minimum  value  A(t)min  = |l-e2(o)J  A(o)  . 

Furthermore  the  time  evolution  of  the  stability. factor  e is 


e2(t)  = 


= + t1 ' '2<o) 
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and  hence  decreases  with  time.  Therefore,  the  stability  of  the  orbits 
is  improved  and  the  mirror  cross  section  of  the  ring  is  reduced  during 
compression. 

IV.  Discussion  and  Conclusions 

The  constrained  or  constant  radius  magnetic  acceleration  of  a non- 
neutral ion  ring  has  two  important  advantages  over  the  free  or  adiabatic 
acceleration.  First,  the  energy  of  the  nonrelativistic  ions  increases 
with  the  square  of  the  magnetic  field  and  second,  the  extraction  of  the 
ions  out  of  the  compression  system  appears  to  be  simpler  than  in  the 
free  compression. 

A magnetic  field  configuration  that  is  suitable  for  constant  radius 

acceleration  of  an  ion  ring  is  the  betatron  field,  that  satisfied  the 

1:2  flux  rule,  i.e.,  the  average  magnetic  field  is  twice  the  magnetic 

field  at  the  orbit  of  the  gyrating  charged  particle.  However,  the 

orbits  of  the  particles  in  a betatron  field  become  unstable  when  the 

self  fields  of  the  ring  are  not  negligible.  Specifically,  the  self 

field  index  ng  £=  qq  (t)fi®(t)J  must  be  less  than  i in  order  for 

the  orbits  to  be  stable.  This  imposses  very  stringent  limitations  on 

the  maximum  current  of  the  ring  as  may  be  seen  from  Eq.  (l8).  For  1 MeV 

protons  forming  a ring  with  an  aspect  ratio  (r  /a)  = 10  and  n = i,  this 

o s 

equation  states  that  the  maximum  current  of  the  ring  cannot  exceed  10  A. 
This  difficulty  is  avoided  by  superimposing  an  azimuthal  magnetic  field 
on  the  betatron  field  (modified  betatron). 

The  Bg  field  improves  considerably  the  stability  of  the  orbits  and 
permits  ng  > > 1,  provided  that  ng  < 2Pq)2  , where  is  the 
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cyclotron  frequency  corresponding  to  the  Bg  field  and  the  cyclotron 
frequency  corresponding  to  the  axial  magnetic  field  at  the  orbit  of  the 
particle. 


In  addition,  the  stability  of  the  orbits  is  improved  with  the  com- 
pression because  the  self  field  index  n^  ^ VQ|  (t)  0Qg(t)  decreases  with 
time  faster  than  the  ratio  (^/g^)2  • The  same  conclusion  is  also 
drawn  from  Eq.  (25).  Since  neither  the  major  nor  the  minor  radius  of 
the  ring  varies  with  increasing  magnetic  field  [See  Eq.  (24)]  and  since 
the  total  number  of  ions  in  the  ring  is  concered,  Eq.  (l8a)  states  that 
ns  * YQ|(t)P  ^g(t).  Substituting  this  expression  for  ng  in  Eq.  (15)  we 
obtain  that  the  field  reversing  factor  C also  decreases  with  increasing 
field  and  is  proportional  to  Y0e(t). 

me  present  theoretical  model  is  based  on  the  assumption  that  the 
radius  of  a typical  gyrating  particle  is  controlled  by  the  external  mag- 
netic field  [see  Eq.  (l)].  Therefore,  it  is  required  that  the  field 
reversal  factor  C <<:  1.  Since  £ decreases  with  time,  the  validity  of 
the  model  is  improved  with  compression.  For  those  applications  requir- 
ing field  reversing  magnetic  field  configurations,  the  field  reversal 
factor  can  be  increased  after  the  acceleration  process  by  adiabatically 
compressing  it  to  a smaller  radius. 

When  n » 1 or  for  arbitrary  n when  n - the  betatron  frequencies 
s s 

ans  u>z  are  equal.  As  a consequence  of  this  fact  the  minor  cross 
section  of  the  ring  retains  its  original  shape,  which  is  assumed  to  be 
circular.  Thus,  the  equivalent  electric  field  in  Eq.  (ll)  is  a good 


approximation,  provided  that  the  aspect  ratio 


H 


>>  1.  As  stated 


above,  for  strong  rings,  i.e.,  when  n >>  1,  stability  of  the  orbits 

s 

requires  that  the  azimuthal  magnetic  field  be  much  greater  than  the 


l4 


f 
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axial  field  at  t = o.  For  this  reason  in  a practical  device  the  proton 
beam  must  be  injected  into  the  system  along  Bg  and  not  along  the  sym- 
metry axis  of  the  system.  It  is  rather  unlikely  that  polarization 
fields  will  introduce  difficulties  when  the  injected  beam  is  space 
charge  neutralized  because  (Bg /Bq)  > > 1. 

The  preceeding  analysis  is  also  based  on  the  assumption  that  the 
ion  ring  is  electrostatically  and  magnetically  noneutral.  For  neutral- 
ized rings  the  effect  of  the  background  electrons  on  the  equilibrium  as 
well  as  on  the  stability  of  the  ring  must  be  included.  The  effect  of 
the  electrons  on  the  ion  ring  depends  on  the  electronic  invironment  of 
the  ring,  that  is,  whether  the  neutralized  ion  ring  is  in  a vacuum  or 
submerged  in  a dense  plasma. 

When  the  ion  ring  is  electrostatically  neutralized  by  background 
electrons  occupying  the  same  volume  as  the  ion  ring,  the  induced 
electric  field  Eg,  will  produce  a zeroth  order  guiding  center  drift  of 
electrons  towards  the  axis.  This  inward  radial  drift  of  electrons  sets 

up  a radial  well  as  a E x B drift  which  tends  to  enhance  the  ion 

~r  ~z 

current.  The  radial  polarization  field,  set  up  by  the  initial  Eg  x 
drift,  acts  to  pull  the  ions  in  radially  if  > > 1.  However,  the 

presence  of  a large  aximuthal  magnetic  field,  such  as  > Bz(t)  re- 
duces the  radial  drift  of  the  electrons  by  the  factor  Bz(t)/Bg  < < 1. 
If,  on  the  other  hand,  the  ion  ring  is  submerged  in  a dense  plasma 
background  of  density  much  greater  than  the  ring  density,  the  radial 
polarization  electric  field  set  up  by  the  Eg  x Bz  drift  of  electrons  is 
shorted  out  by  the  plasma  electrons.  Hence,  the  ion  ring  maintains  a 
fixed  average  radius  as  the  axial  magnetic  field  is  increased. 
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The  equilibrium  and  stability  of  ion  rings  in  the  presence  of  a 
space  charge  neutralizing  electron  background  is  under  investigation 
with  a particle  simulation  computer  code. 

Appendix  A 

In  this  appendix  a general  discussion  of  the  solution  to  the  orbit 

equations  given  in  Eqs.  (l4a,b)  is  presented.  Since  V 0 , 2 and  u>  2 

oo  r z 

are  assumed  to  be  weak  functions  of  time  the  solutions  to  (l4a,b)  can 
be  written  in  the  form 

• A 

6r  = a e'1  J *(t')dt'  , (Al) 

o 

6z  = b e"1  J1  ^(tOdt'  , (A2 ) 

o 

where  a,  b are  constant  and  u>(t*)  is  in  general  complex.  Substituting 
(Al)  and  (A2)  into  Eqs.  (l4a,b)  and  noting  that  |d>!  « |<*>2|  we  obtain 


where  u>2  = (no/Yo0)2  (l-n-ns)  and  u>2  = (^o/yo6)2  (n-ng). 

Combining  (AJ ) and  (A4)  and  keeping  only  the  lowest  order  terms  in 
the  time  derivatives  of  u>  and  Yog  the  following  dispersion  relationship 
correct  to  first  order  in  u>/u>2  and  y/y2  is  obtained 


1 6 


I 


((JO2  - U>2)  («»2  - CO2)  - 102 


(-) 

\\e) 


+iu>  2iu2  - co2  - u)2  ( — + = 

r ZJ  V°  Yo0/ 


(A5) 


We  now  let  uo  - , 


where  < < jco^^is  such  that  it  satisfied  the  following 

zeroth  order  time  dependent  dispension  relation 


(A6) 


(u>(°))2  - <o2J  (u/°^)2  - U)2  - (u'(°))2  (ni/vo0)2  = o . (A7) 


Substituting  A(6)  into  (A5)  and  using  (A7)  we  arrive  at  the  following 
expression  for 


2(u>(o))g  - v2  + *(o)/u>(o)+  Yoe/v0e]  ? (a8) 


2(ur^)2  - v2 


where  v2  = u>2  + u>z2  + M^i/2YQg)2  and  is  given  by  the  solution 

of  (A7).  The  four  solutions  of  Eq.  (A7)  are 


to 


i?2  - + (v  | + ^ v4  - lm*  / ^T" 


<A9) 


k I 
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and 

= ± (v2  ~ 

V ^ 

The  relations  in  (A8-10)  together  with  (Al)  and  (A2)  formally  describe 

the  linear  evolution  of  the  perturbed  orbits. 

Noting  the  definitions  of  u>  2 and  u»  2,  we  see  from  Eqs.  (A9)  and 

I*  z 

(A10)  that  stability  is  achieved  if  the  following  inequalities  are 
satisfied 

ns  < i |"l  + (ni/k0)2]  ’ (All) 

and 

(l  - 2ns)  + (^i/^0)2  1 > 4 (n-ns)  (l-n-ns)>  0 • (A12) 


{' 


- kw  ^ : 
r z 


1.  . 


(A  10) 


In  the  special  case  where  n^_  — —o  stability  is  implied  if 

o < n < 1 . (A13 ) 

The  inequality  in  (A13)  is  the  well  known  stability  condition  for  a 
single  particle  in  the  absence  of  a Bg  field.  If  = o and  n ^ o, 
the  orbits  are  stable  if 

ns  < § , (Al4) 

is  satisfied  together  with 

n(l-n)  < i (A15) 
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However,  when  = o and  ^1  ^ o the  condition  for  stability  is  the  same 

as  when  n = ^1  = o,  that  is,  the  external  field  index  must  satisfy  (A13), 
s 

Hence,  in  the  absence  of  self  fields  the  regime  of  stability  is  not  in- 
creased by  the  introduction  of  an  azimuthal  magnetic  field.  A situation 
of  somewhat  more  practical  importance  is  when  the  self  field  index  is 
much  greater  than  unity,  i.e.,  no  > > 1 and  the  external  field  index  is 
of  the  order  of  unity,  |n|  = 0(l).  For  this  case  the  stability  condi- 


tion is  given  by 


n < k 


(A16) 


When  n^  >>  1,  the  radial  and  axial  betatron  frequencies  are  approxi- 
mately equal,  i.e.,  u>2_  tw2=  u)2  . The  stability  condition  in  (Al6) 

can  then  be  stated  in  the  form  e2  < 1 where  e2  = ^/Pi/2Y  j )2 . To 
obtain  explicit  expressions  for  and  from  Eqs.  (a8)  and  (A9)j 

we  consider  the  case  where  n > > 1,  |n|  _ 0(l)  and  e2  < < 1,  which 

s 

corresponds  to  a dense,  stable  ring.  The  zeroth  order  frequencies  in 
(A9)  and  (A10)  for  this  case  are 


(2  + e2/2 


(A  17) 


3 4 


(A18) 
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Substitution  of  (A17)  and  (Al8)  yields  the  single  first  order  frequency 


Fig.  1 — Schematic  drawing  of  an  apparatus  that  can  be  used  for  the  magnetic  acceleration 
of  a proton  ring.  The  Bfl  field  may  be  produced  by  the  current  IQ  flowing  in  the  center  con- 
ductor is  independent  of  time.  In  the  laboratory,  the  B0  field  may  be  produced  convenient!'' 
with  toroidal  windings.  The  Bz  and  Br  fields  vary  with  time  and  are  produced  by  the  two  sets 
of  concentric  coils.  In  the  proposed  scheme,  not  only  the  major  radius  rQ  but  to  a good  ap- 
proximation also  the  minor  radius  of  the  ring  remains  constant  during  acceleration. 
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